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The general features of the evolution of the Yukawa coupling constants and seesaw operators in 
the universal seesaw model with detMi? = are investigated. In the model, not only the magnitude 
of the Yukawa couphng constant (Y^ )33 in the up-quark sector but also that of (Y^ )33 in the down- 
quark sector is of the order of one, i.e., (Y2')33 ~ (Y/)33 ~ 1. The requirement that the model 
should be calculable perturbatively, i.e., |^jP/47r < 1, puts some constraints on the values of the 
intermediate mass scales and tan/3 (in the SUSY model). 



I. INTRODUCTION 

Recently, there has been considerable interest in the 
evolution (energy-scale dependency) of the Yukawa cou- 
pling constants of quarks and leptons. If we intend to 
build a model which gives a unified description of quark 
and lepton mass matrices, we cannot avoid investigat- 
ing evolutions of the Yukawa coupling constants. The 
recent study on the quark masses and mixings has been 
given, for example, in Ref. Especially, recently, the 
evolution of the neutrino seesaw mass matrix has been re- 
ceived considerable attention (for example, see Ref. ||]) 
in connection with the energy-scale dependence of the 
large mixing angle. 

As one of such unified models, there is a non-standard 
model, the so-called "universal seesaw model" (USM) [|. 
The model describes not only the neutrino mass matrix 
Ml, but also the quark mass matrices M„ and Md and 
charged lepton mass matrix by the seesaw-type ma- 
trices universally: The model has hypothetical fermions 
Fi {F = U, D, N,E; i = 1, 2, 3) in addition to the conven- 
tional quarks and leptons fi (/ = u,d,v,e; i = 1,2,3), 
and these fermions are assigned to /l = (2,1), fn, = 
(1,2), Fl = (1,1) and Fn = (1,1) of SU(2)lx SU(2)«. 
The 6x6 mass matrix which is sandwiched between the 
fields (Jl.Fl) and ifR,FR) is given by 

^6X6 ^ M ™^ (1.1) 

where and tur are universal for all fermion sec- 
tors (/ = u, d, V, e) and only Mp have structures de- 
pendent on the flavors F. For Al < A^ ^ A5, where 
Al = O(mL), Kr = 0{mR) and A5 = 0{Mf), the 3 x 3 
mass matrices Mf for the fermions / are given by the 
well-known seesaw expression 

Mf ~ -rriLMp^m^j^ . (1.2) 

Thus, the model answers the question why the masses 
of quarks (except for top quark) and charged leptons are 



so small compared with the electroweak scale A^, (^ 10^ 
GeV). 

Recently, in order to understand the observed fact 
rrit Ai (mj is the top quark mass), the authors have 
proposed a universal seesaw mass matrix model with an 
ansatz detMp = for the up-quark sector {F = U). 
In the model, one of the fermion masses m[Ui) is zero 
[say, m{U^) = 0], so that the seesaw mechanism does not 
work for the third family, i.e., the fermions {uj,l,U^r) 
and (J73l,U3_r) acquire masses of 0{mL) and 0{niR), 
respectively. We identify [u^LiU^j^) as the top quark 
{tL,tR). Thus, we can understand the question why only 
the top quark has a mass of the order of A^. 

Our interest is as follows: In the conventional model, 
the Yukawa coupling constants yf of the fermions / 
are given by ?// = mf/{(l>1). Only the Yukawa cou- 
pling constants yt of the top quark t takes a large value 
yt = ^t/{4>^L) ~ 1- The other Yukawa coupling con- 
stants yf are sufficiently smaller than one. On the con- 
trast to the conventional model, in this USM, the ma- 
trices = Y[{(jPj^) are universal for all fermion sectors 
f = u,d, e, V, i.e., = Y/ = Yf = Y^. Therefore, when 
(Y2')33 is of the order of one, the other (y/)33 will also be 
of the order of one. We are afraid that in such a model the 
Yukawa coupling constants have Landau poles at energy 
scales lower than a unification energy scale ^ = Kx (so 
that the model causes "burst" of Yukawa coupling con- 
stants before going to the unification energy scale). One 
of our interests is to investigate whether such a model 
can provide or not a set of reasonable parameter values 
under the conditions that the Yukawa coupling constants 
(and also the seesaw operators m^Af^^?/!,^) do not have 
the Landau poles below /i — Kx- 

We also take an interest in the "democratic" USM , 
which is an extended version of USM and has success- 
fully given the quark masses and the Cabibbo-Kobayashi- 
Maskawa (CKM) Q matrix parameters in terms of the 
charged lepton masses. However, the study is only phe- 
nomenology at the energy scale = mz (jnz is the neu- 
tral weak boson mass). Since the model is one of the 
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promising models of the unified description of the quark 
and lepton mass matrices, it is important to investigate 
the evolutions of the mass matrices in the USM. 

The democratic USM is as follows: 
(i) The mass matrices rriL and niR have the same struc- 
ture except for their phase factors 



(1.3) 



where k is a constant and are given by 

= diag (^zi exp(i(5(), Z2 exp(i(52), exp{i6l)^ , (1.4) 

with zf + Z2 + z'^ — 1. (The fermion masses m( are in- 
dependent of the parameters S{ . Only the values of the 
CKM matrix parameters \Vij\ depend on the parameters 



(ii) In the basis on which the matrices m£ and m}^ are 
diagonal, the mass matrices Alp are given by the form 



moA(l + 3&/X), 
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(1.5) 



(1.6) 



(iii) The parameter bf for the charged lepton sector is 
given by be = 0, so that in the limit of k/A <C 1, the 
parameters Zi are given by 



Z2 



Z3 



1 



^/rr^e 



Tot 



(1.7) 



Then, the up- and down-quark masses are successfully 



given by the choice of b^ 



-1/3 and bd 



{(3d = 18°), respectively. The CKM matrix is also suc- 
cessfully obtained by taking 



"3 



(1.8) 



However, when we take the evolution of the Yukawa 
coupling constants into the consideration, we should con- 
sider that the assumptions (i) and (ii) are required not at 
the electroweak energy scale /i — , but at a unification 
energy scale ^ = Ax, i.e., the assumptions (i) and (ii) 
should be replaced with 



Yl{Ax)^Yl{Ax)^^ij,Zf 



and 



Y^{Ax)^esi^ + SbfX) 



- cS 



(1.10) 



(1.11) 



respectively, where mass matrices to^, to^ and Mp are 
expressed by 



Y 



y/($°) 



respectively, and 



(1.12) 



(1.13) 



and 0L, (1)11 and $ are Higgs scalars whose vacuum ex- 
pectation values (VEV) break SU(2)l, SU(2)fl, and an 
additional U(l) symmetry U(l)xj respectively. (For sim- 
plicity, we have assumed that the values of (02,); ^'^d 
(<f>°) are real.) 

Another interest in the present paper is to check 
whether or not the phenomenological study in the previ- 
ous paper Q| is still approximately valid under the evolu- 
tion of the Yukawa coupling constants. For example, the 
model with be = and 5„ = —1/3 has led to the relation 



rrie 



3 TOe 

4 TO,, 



(1.14) 



almost independently of the value of the seesaw suppres- 
sion factor At/A. One of the reasons to taking the value 
of 6/ in the up-quark sector as 6„ = —1/3 exists in the 
successful relation (1.14). Therefore, we have interest 
whether the relation (1.14) still holds even when we take 
the evolution into consideration. 

Besides, even apart from such phenomenological inter- 
ests, it is very important to investigate the general fea- 
tures of the evolution of the Yukawa coupling constants 
in the universal seesaw model with detMp = 0, because 
in the present model one of the fermions Fi does not de- 
couple from the theory at /i < Ag, so that the evolution 
shows peculiar behavior in contrast with the conventional 
seesaw model. 

A similar study has been done in Ref. ^ by one of the 
authors (Y.K.). However, in Ref. ||], instead of the see- 
saw operators which will be defined later in Eqs. (3.8) 



corresponding to to^M^ 



the evolution of the seesaw 



forms of the Yukawa coupling constants Yl{Yg) ^{Y^Y 
were investigated by calculating the Yukawa coupling 
constants y/ , Y^ and Yg individually under the assump- 
tion that the heavy particles with the masses of the order 
of As = ($°) do not contribute to the evolution of Y^ 
{A — L,R,S) below /i = Ag. In the present paper, 
we will calculate the evolution of the Yukawa coupling 
(1-9) constants y{ above /i = A5 and that of the seesaw op- 



erators K-^ below /i 
in Sec. 



Ill 



A5, except for (y/)i3 as discussed 



In Sec. H we will discuss an additional symmetry 
which is introduced for the purpose of preventing that 
the fermions F acquire the masses Mp at the energy 
scale /i ~ As- In Sec. Ill, we will give the general formu- 



lation of the evolution of the seesaw mass matrices with 
detMjj = 0. In Sec. IV, we give the explicit coefficients 
of the renormalization group equations. In Sec. 0, we 
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discuss the evolution of an extended version of the USM, 
the "democratic seesaw model" [HJ^. The numerical re- 
sults for a non-SUSY model and for a minimal SUSY 

respectively. 



VII 



It 



model are given in Sees. VI and 
will be emphasized that the energy scale dependencies 
in the SUSY model are quite different from those in the 
non-SUSY model. The evolution of the neutrino mass 



matrix is given in Sec. VIII . It will be showed that, dif- 
ferently from the conventional seesaw model, the present 
neutrino mass matrix is form-invariant below /i = A5. 
Finally, Sec. IX will be devoted to the conclusions and 
remarks. 



II. U(l)x SYMMETRY 



In the present model, the gauge symmetries are broken 
as follows: 



SU(2)^ X SU(2)^ X U(l)^^ X SU(3), x U(l)^ 
i M = A5 
SU(2)^ X SU(2)^ X U(l)^« X SU(3), 

i ^i = ^R (2.1) 
SU(2)^ X U(l)^ X SU(3)^ 
i /i = Al 
U(l)_ X SU(3), . 

Here, the symmetry \]{l)x, which is spontaneously bro- 
ken at the energy scale /i = A5, has been introduced for 
the purpose of preventing that the fermions F acquire 
the masses Mp at > A5. Hereafter, we call the ranges 
< IJ' < Afl, Afl < fi < As, and A5 < /i < Ax as the 
ranges I, II, and HI, respectively. In the present paper, 
the energy scale Ax does not always mean a gauge unifi- 
cation energy scale. We assume that at the energy scale 
Ax the mass matrices (Yukawa coupling constants) take 
simple forms discussed in the previous section. 

The Yukawa coupling constants Y^, and Yg are 
defined as follows: 



+ Ymj^B.iRULj + Y^^^qj^^^rDlj + Y^^~iRi4>RNLj + Y^^/r,c^rElj 

+ Y^^Pm'^Ur, + Yi^^Du^^DR, + Y,^,^Nu<i>NRj + Y|^^^L^'^^ Er, + h.c. 



(2.2) 



where 

iL/R 
(kh/R 



L/R 



L/R. 



tLIR^ 



^LjR 



L/R. 



(2.3) 



L/R. 



From Eq. (2.2), the U(l)x charge assignment should 
satisfy the following relations 



X{Ur) = X{Ul) - X($) , 
X{Dr)=X{Dl)+X{<^) , 



X{qL) = -[X{UR)+X{DR)] 
X{qR) = \[X{UL)+X{DL)] 



(2.4) 



(2.5) 



X{<j>L) 

X{cj>R) 



[X{Ur) - X{Dr)] 
[X(Ul)-X{Dl)] 



(2.6) 



for quark sectors, and equations similar to Eqs. (2.4) - 
(2.6) for lepton sectors f = v,e. For simplicity, in the 
present paper, we choose 



X{qL/R) - X{eL/R) - , X($) = +1 



(2.7) 



Then, the quantum numbers of the fermions / and F 
and Higgs scalars 0l, 4>r and $ for SU(2)^ x SU(2)^ x 
U(l)^^ X U(l)j(- are given in Table | 

Note that the quantum number of the fermion is 



identical with that of the fermion TV' 



R 



{NrY = CNf 



Therefore, the neutral fermions Nl and Nr can acquire 
the following Majorana mass terms at /i = A5: 



H 



Majorana 



Ys%Nl^N2^+Ys%Nr,i.r, 



h.c. 
(2., 



Then, the neutrino mass matrix is given as follows 



Vl V^r Nl Nj^ 
















m 


T 

mi 






Ml 
Ml 



rriL \ 


Md 

Mr) 



( 

VR 

Nl 
\Nr 



(2.9) 



Y/($) 



where Md = Y^{<^), Ml = Y^{^) and Mr 
Since 0{Md) ~ 0{Ml) ~ 0{Mr) » O(m^j) > O(mL), 
we obtain the mass matrix M^, for the active neutrinos 

VL 



-mLM„^mT 



(2.10) 
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III. GENERAL FEATURES OF THE 
EVOLUTIONS 

In the present section, we give a general formulation 
of the evolution of the seesaw matrix with detMp = 0. 
The evolution of the neutrino seesaw mass matrix is well 
known. However, in such a model with detAf^? = as the 
present model (the democratic seesaw model), a careful 
treatment is required. 

Without losing the generality, we can express the 
Yukawa coupling constants Y/ and (/ = u, d, v, e) 
as 

= , = , (3.1) 

where Z^{^) {A — L,R) are defined by 

Z{{fj,) = di£Lg{z{^{n),z{2{fi),z{^{^i)) , (3.2) 



+ 1^2(^)1' + 1^3(^)1' 



(3.3) 



on the basis on which (fi) are diagonal. In the present 
model, the word "universal" means the following initial 
conditions 



-Li 



(Ax)| 



= l4(Ax)| 



(3.4) 



(3.5) 



for all fermion sectors f = u,d,v^e universally. 

In the range III (A5 < fi < Ax), the evolutions of the 
Yukawa coupling constants , and Yg are given by 
the one loop renomalization group equations (RGE) as 
follows: 

16n^^^(T{~G{ + H{)Y{ , (A = L,i?,5), 

(3.6) 

where t = log/i, and t{, G;^ and h{ {A = L, R, S) de- 
note contributions from fermion loop corrections, vertex 
corrections due to the gauge bosons, and vertex correc- 
tions due to the Higgs boson, respectively. Note that the 
matrices and G]^ are proportional to the unit ma- 
trix. As stated in the next section, the coefficients 
{A = L, R, S) take diagonal forms on the basis on which 
Y^ are diagonal. Therefore, if we take a basis on which 
y/ (and Y^f ) or Yg are diagonal at /i = A^, then the 
Yukawa coupling constants y/ (and Y^) or Yg can keep 
the forms diagonal in the range III. Sometimes, the basis 
on which y/ (and Y^) are diagonal is useful, but some- 
times, another basis on which Yg are diagonal is useful, 
as we discuss later. 

In the present model, it is assumed that we can choose 
a flavor basis on which y/(Ax) are simultaneously diag- 
onal for all f = u,d, ly, e. Then, on this basis, since the 



Yukawa coupling constants Yg {fj) can keep the forms di- 
agonal in the range III, we can find that all Yg are diag- 
onal 'dt fi — As- We can denote those as 



y/(As) = diag{y(g,ylg,y(g) 



(3.7) 



At the energy scale /i = As, the fermions Fi (except 
for C/3) acquire the heavy masses {Mp)ii = y{g{^'^)- In 
the conventional seesaw model with det Mp ^ 0, the en- 
ergy scale behaviors of the fermion masses in /i < A5 are 
described by evolutions of the following operators 



Y^l-iY[)AY!,)% 

k=l VkS 



(3.8) 



and 



k=i 



jk ■ 

(3.9) 



(Hereafter, for convenience, we will denote the Yukawa 
coupling constants Yg^ in the Majorana mass matrix 
Mn — Y(?(<I>°) as Yg.) The quark and lepton mass ma- 
trices Mf are given by 

M/ = X/(0°)(0Oj)/(<l>°) , if = u,d,e), (3.10) 



(3.11) 



As explicitly shown in Sec. iV, the evolutions of the op- 
erators Rf are described by the one- loop RGB's with the 
following forms 



167r^ 



,d_K^ 
IT 



167r^ 



.dK" 
IT 



{f^u,d,e), (3.12) 



= (Tj^ - G'i,) K-^ + H^K^K-^ + K'^H^/^ , 



(3.13) 



where T^, G;^ and {h{^j^, H^.^^) denote contributions 
from fermion loop corrections, vertex corrections due to 
the gauge bosons, and vertex corrections due to the Higgs 
bosons (t)L and (pn, respectively. 

However, in the seesaw mass matrix with detM^? = 0, 
since one of the eigenvalues of Yg (/ = u) is zero (say, 
Vss = 0), we must calculate the following operator 



(K-),,^ Yl{YS)-^Y, 



"t 



T.-r{yL)^u{YR)% 

fc=l ^kS 



(3.14) 
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di&g{y'^g,y^s). Note that the matrices Yu rrif = {(!>%) J \{Y^)i3? + |(ifl)23p + |(^fl)33P • (3.17) 



where Yg 

and 17 (^r) in Eq. (3.14) are 2 x 2 and 3 x 2 matrices, 
respectively. 

Note that in (3.14) we have taken the sum over k = 1 
and 2 only. In the range II, the evolutions of the Yukawa 
coupling constants Y^is ™d ^mz (* = 1:2,3) are still 
described by the equation (3.6). At the energy scale 
fjL = Kr, we obtain a new mass term 

Hmass = Y.^Yii)*^UL3Um{^l) . (3.15) 



By defining a mixing state 

^/\{YRU? + \{Y^U? + \{YRU? 



^R3 



(3.16) 



we obtain a mass me of the fourth up-quark t' = 



Similarly, in the approximation in which the terms sup- 
pressed by y^g and y2s are neglected, the mass mt of the 
third up-quark (i.e., top quark) t = {tL,tR) = (u^s, Urs) 
is given by 



mt ^ (<^i)^|(l?)l3P + |(1?)23P + |(1?)33P , (3.18) 

where 

{Yi:)hULl + {Yi:)hUL2 + (Y^hULS 



Ul3 - 



^/\(YE)WT\mwT\(Y^^ 



. (3.19) 



More precisely speaking, the masses {mu,mc,mt,mt') 
are obtained by diagonalizing the following mass matrix 



-{k/X)K, 

'*-'^iJ13 



11 


-(«/A)xr2 




u 
21 






'U 

31 










'*-'^fl33 



'-L23 
-'L33 



(3.20) 



which is sandwiched by the fields {uli,ul2,ul3,Ul3) 
and {uri,ur2,ur3,Ur3), where k = {(I)r)/{<I)l) and 
A = {<i>°)/{(t)l) as defined in Eq. (1.13). Of the Yukawa 
coupling constants {Y^)ij and {Y^)ij^ the twelve com- 
ponents iYj^)ik and {Y^)ik {i = l,2,3;j = 1,2) are ab- 
sorbed into the operator defined by (3.14), while the 
rest (YL)ii and {Y^)i'i are still described by the equation 
(3.6). 

Finally, we denote the effective Hamiltonian in each 

range: The effective Hamiltonian Hf^H in the range III 
(Ajf > /X > A5) is still given by the form (2.2), and 
in the range II (As > /x > A^) and Hl^^ in the range I 
(Afl > IJ,> Al) are given by 



H^t=J2YMQLi<l>LUR3) 

i=l 
3 

+ Y.'^Ri3iQRi^RUL3) 
i=3 



and 



HL = Y.^Li3{UiHURz) 
i=l 

+ E ^^^m^u4>LUR,) 



+ E 



''''^ KlMLicl>LeRj) + h.c. 



(cf>0) 



respectively. 



(3.21) 



(3.22) 



IV. COEFFICIENTS OF THE RGE 

In the present section, we give the coefficients of the 
renomalization group equations (RGE) (3.6), (3.12) and 
(3.13). 
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A. Evolution in the range III 

In the non- SUSY model, the terms t{, G{ and h{ 
{A = L, R, S) are given as follows: 

rpu fj^d rpu rpe 

= 3TV (yXYX^ + YiYi^) + Tr (yXYX^ + YXY^^) , 

(4.1) 



167rVln 
dt 



[T[-Gi + Hl)u-{Ti,-Gi + Hih 



(4i 



The case g2L = .92/? is likely in the L-R symmetric 
model. For convenience, in the numerical evaluation in 
the present paper, we will take 521, (m) = 92r{^i) in the 
range III (As < < A^). 



5 

s' 

9 

s' 

45 



Ga = ^9l + \91a + ^9% 



HI = -Hi = l {YXYf - YiYi^) , 



^A — -^A 



^^^YXYX^ 



YaYa 



where A = -L, i?, and 



(4.2) 



3Tr (Fs"F^t ^ ydy^^t j ^ ^r {y^Y^'^ + YIY'g^ 



(4.4) 



B. Evolution in the ranges I and II 

In the ranges I and II, all the fermions Fi except for 
are decoupled from the equation (3.6). In the present 
section, we will take the diagonal basis of Mp- Therefore, 
it is convenient that we define a spurion 




5= ( 
1 



(4.9) 



(4.3) Then, the surviving Yukawa coupling constants {YX)iz 
are expressed as {YXS)ij = {YX)iz5^j. The evolution of 
YXS is still described by the RGE (3.6) by substituting 
YXS for . Here, the terms Tjf, G\ and H^{A^L,R) 
are expressed as follows [Y^ {f = d,e,v) are already ab- 
sorbed into the operators K^]: 



TX = 3Tr {yXSYX^^ 



(4.10) 



G^s 
G% 
G-s 
G| 



4^1 + Sgl + -g\ 
9l + ^9l + \9\ , 



+^9x, 
Q9! 



Iff ■w-f-w-f'f 



:9x 



if = d, ly, e) 



(4.5) 



(4.6) 



The coefficients T^, G{ and H{ in the minimal SUSY 
model are given in the Appendix. 

As seen from Eq. (4.6), since the matrix is diagonal 
on the diagonal basis of Mp{Ax), the Yukawa coupling 
constants Yg (fj.) can keep the forms diagonal. Similarly, 
when we choose the diagonal basis of Mj:,(Ax) [Mji{Ax)], 
the matrices Y[(ji) [Y^(/i)] keep their forms diagonal. 



For a model with g2L(A*) = 92r{iJ') and Y^ (/x) 
at n = Ax , we can assert that 



(4.7) 



in the range III (Ag < ji < Ax), because on the diagonal 
basis of Yl we obtain 



G\ = + \91a 



HI = -YXSY^^ 



A > 



{A = L,R) in. the range II, and 

T^ = 3Tr(Y^SY^''^ , 



Gl = '^gl + I9L + 89I , 



HI = \y}:SY-^ , 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



in the range I. Here, the coupling constant gi = giLR in 
the range II is that for the U(l) operator {1/2)Ylji which 
is defined by the relation 



(4.16) 
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for the symmetry SU(2)^ x SU(2)^ x U(l)^^, while the 
coupUng constant gi = giy in the range I is that for the 
U(l) operator {1/2)Y which is defined by the relation 



(4.29) 



3 ' 2 ' 



(4.17) 



for the symmetry S\J{2)j^ x U(l)y, and they are con- 
nected by 

aem(Ai) = a^li^L) + ^a-2«(Ai) , (4.18) 

^a-^iAn) = a-^{AR) + ^a-ljiiAR) , (4.19) 

where ai = gf/Air. 

Similarly, the terms T^, G;^ ^kl ^^'^ ^kr if ~ 
u, d, e) are given by 

= = = 3Ti (y^ S Y^^ +Y^SY^^^ , (4.20) 



= Gi = Ig! + lal, + \gl^ + Sgj , 

9 9 9 9 9 9 

Gk = + ^dk + -^qIr , 



(4.21) 



Ili;,A = H'KA = ^yXSYX^ , if^^ = 0, {A = L,R), 

(4.22) 

in the range II, and 

T^ = T^=T^ = 3Tr (y^SY^^'^ , (4.23) 



(4.24) 



GK=^^9l + l9L+8gl, 

G'k = ^9f + l9L + 89l, 

_ 45 2 82 
Gr - -^91 + -^92L > 



HlL=H'ki^ = lY]iSY-\ Hi^ = 0, f = u,d, 



(4.25) 
(4.26) 



in the range I. 

The terms Tj^, and H^j^ have rather simple forms 
in contrast with those in the conventional neutrino see- 
saw model, bccaiise the partners of the fermions which 
couple to the Higgs scalar are not Jr, but Fr which 
are already decoupled at /z = Ar: 



Tj^ = 6Tr (y^'SY^t) , (4.27) 
= 3gk , (4.28) 



in the ranges I and II, where \hl is the coupling constant 
of the Higgs scalar (pj, defined by 



H4, = ^XHL{(t)L(f>Lf , 



(4.30) 



and the mass of the physical Higgs scalar is given by 



(4.31) 



The similar coefficients in the minimal SUSY model 
are given in the Appendix. 



V. CASE OF THE DEMOCRATIC SEESAW 
MODEL 

In the democratic seesaw model, on the diagonal basis 
of Y[{Ax) and Y^(Ajf ), the Yukawa coupling constants 
of heavy fermions Yg{Ax) are given by the democratic 
form (1.11). Since on this basis the Yukawa coupling 
constants Yg keep the forms democratic: 



^/(m)=C|(m)(1 + 36/(m)X) 



(5.1) 



we will call this basis the "democratic basis of Mj?" here- 
after. On the other hand, if we take a basis on which Yg 
are diagonal, i.e., the matrix forms are given by 



Yi{p)=es{i^){i+ibf{n)x) 

X = AXA^ = diag(0, 0, 1) , 



A = 



1 



V3 V3 




(5.2) 



(5.3) 



(5.4) 



Especially, on this basis, the Yukawa coupling constants 
{Yg)ii and {Yg)ii of the fermions Ei and Ui satisfy the 
relations 

[y|(M)]ii = [F|(/x)]22 = [F|(m)]33 = ^I(m) , (5.5) 



[>s (m)]ii = [Y^{t^)h2 = Csif^) , [Y^{f^)hs = , (5.6) 

in the range III (A5 < fi < Ax), i.e., 

Will) = , huifj) = -1/3 . (5.7) 

On the other hand, on this basis, the Yukawa coupling 
constants Y^ip) and yJ^{ii) are not diagonal. However, 
we can easily obtain their diagonal forms by A^yI{ii)A 
and A'^yI{ii)A. 

At the energy scale /i = A5, the fermions Fi (except 
for U3) acquire the heavy masses {Mp)ii. Therefore, for 
11 < As, the operators are given as follows: 
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fe=l,2 



(5. 



(5.9) 



{K^).,= YliYir'Y^^ 



{YtUiY^Tjk 



2,3 



(5.10) 



-J CS(As) 



, fc=l,2 



1 



1 + 36, (As) 



(5.11) 



In Eq. (5.11), we have assumed that the structure of 
the Majorana mass term MYi{l^s) = ^/(As)($°) for the 
neutral fermions has a structure similar to the Dirac 
mass matrices Mp which is given by Eq. (5.1). 

Since the Yukawa couphng constants Yj^{^) {A = L, R) 
in the range III keep their forms diagonal on the demo- 
cratic basis of Mp, it is convenient to express ^4 (a*) as 
follows, 



(5.12) 



where the diagonal matrix Z^{ii) is given by Eq. (3.2). 
Then, the matrix Y^ on the diagonal basis of Mp is given 
by 



where 



Yl{^^) = ^{{^^)zf{^^) 



(5.13) 
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3(Z2 + Zi) -\/3(z2 - Zi) 

-V3{Z2 - Zl) 4Z3 + Z2+ Zl 

-VG{z2 - Zl) -\/2(2z3 - Z2 - Zl) 



-\/6(z2 - Zl) 

-V2{2z3 -Z2- Zl) 

2(Z3 + Z2 + Zl) 



(5.14) 



(we have dropped the indices A and /, and for simplicity, 
we have taken 6{ = 0). Although the Yukawa coupling 
constants Y^ and Y^ in the range II and F/ in the range 
I have the physical meaning only for the one column ma- 
trix components (l^)i3 (« = 1,2,3), we still use the ex- 



pressions (5.12) and (5.13), because the matrix K'^{ijl) 
(Al < ^ < As) which is proportional to Yl{ij)Y^ {^) is 
still diagonal on the democratic basis of Mp as discussed 
in Sec. IV, so that we regard that Y'^{^) is also "diago- 
nal". Then, the top quark mass mt{fi) is approximately 
expressed as 



mt(^) 




Li\ 



(5.15) 



The expression (5.15) is valid in the whole ranges Al < 

M < Ax- 
Since 



™(u)+rn (.)^2CE(M)MM<<M) 



(5.17) 



3 2 



i=l k=l 

we obtain 



from Eq. (5.8). Note that the expression (5.17) is valid 
only in the range III. In the ranges I and II, the ra- 
tio Cl^rI^S behaves as a operator K'^{n) which obeys 
Eq. (3.12). From Eq. (5.17), the ratio rric/mt is given by 
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(5.18) 



Since Hf^j^ = H^^j^ = in the ranges I and II, the form 
of K'^(fj,) is invariant in the ranges, i.e., 

Z£(Ai)Z]^t(A^) ^ ziiAs)Z^JiAs) , (5.19) 

especially, since 

Ziip) = Z^(/i) ^ Zfi^,) , (5.20) 

for a model with .g2fl(A7?) = 52L(Afl), we obtain 

Z-(Al) = Z'{As) . (5.21) 

Therefore, in preliminary evaluations prior to fixing the 
final values of the parameters, we will sometimes use the 
values of Zi{mz) which are obtained from the observed 
charge lepton masses mf{mz) by using Eq. (1.7) instead 
of the values of Zi{Ax) which are defined in Eq. (1.9) as 
the initial condition at /i = Ax- 



VI. NUMERICAL RESULTS IN THE NON-SUSY 
MODEL 



We define 



Ar = 



Ar = 



'-'R.I ' 



As = (<f °) 



(6.1) 



However, for convenience, in the numerical evaluations, 
instead of physical quantities at /i — A^, we will use those 
at /X = mz (mz is the neutral weak boson mass). 

First, in order to overlook the behavior of the Yukawa 
couphng constant Y/(/Lt), we illustrate the behavior of 
^£(/^) in the non- SUSY model in Fig. |. Here, we have 
used the approximate relation (5.15) and the input values 
mtimz) = 181 GeV and {(j)^) = 174 GeV: 



mtijnz) 



= 1.80. 



(6.2) 



In other words, the behavior of £,^(p) corresponds to that 
of mt(^) because of = {mtin) /■mt{mz))£.l{mz). 

In the ranges I and II, since the terms and are 
expressed only in terms of Y^^SYj^^ the evolution of the 

factor |CE(m)P = iTr[Y^SY^^] is described by the equa- 
tion 



u\2 



= 2 



Gl 



ICE I 



(6.3) 



However, in the range III, the terms Tf^ and H'f^ contain 
other factors Y[y[^ in addition to Y^Y^\ so that the 
evolution of ^£ cannot be expressed so simply such as 
(6.3). For the evaluation of in the range HI, we have 
tentatively substituted the values Zilmz) given by (1.7) 



for the initial values ^^(Ax). For simplicity, as we dis- 
cussed in (4.7), we have taken as 52l(A_r) = g2i?(A_R). In 
Fig. |l|, the ratio As/Ar has been taken as As/ Ar = 107, 
which has determined from the fitting of the observed 
ratio rat/rric as we discuss later. The behavior of £r(/^) 
is insensitive to the ratio As/Ar. As seen in Fig. |l|, in 
a case with a lower As [As < 10^ GeV), CtiiA has the 
Landau pole below /i = Ax, so that the case is ruled out. 
On the other hand, a case with a higher As [As > 10^^ 
GeV) causes ai(/i) — > cx) at — > A5, so that the case is 
also ruled out. 

Taking account of the behavior of C£(/i) shown in 
Fig. 1^, as a trial, we take 



Ax = 2 X IQis GeV 



(6.4) 



which is known as the unification energy scale in the min- 
imal SUSY model. (However, in the present paper, we 
do not consider the gauge unification.) As a value of A5, 
we tentatively take 



As = 3 X 10^3 GeV , 



(6.5) 



which leads to the mass-squared difference Amgj = 
TO^g — mj^ ~ (10~^ — 10~^)eV^ as we demonstrate later. 
For the values (6.4) and (6.5), we obtain £,l{Ax) = 1.2. 

Next, we determine the values of Cs(^) ^s/-^R- 
Since we have already obtained the value ^£(Ax) = 1.2, 
it seems that we can fix the value of £,g{A)As / Ar from 
the observed value of nit {'mz)/mc{mz) because of the re- 
lation (5.18). However, the value of QiAx) [also ^^(Ax)] 
is sensitive to the value of C5(As) [Cs(^'S')] i^^ other 
words, a small deviation of £,g{As) causes a large de- 
viation of ^^(Ax)]. Therefore, we cannot fix the values 

£,s{Ax) unless we put a tentative model for and 
The basic assumption in the universal seesaw model is to 
consider that the mass matrices tol and tur in Eq. (1.1) 
are "universal" (common) for all fermion sectors (quarks 
and leptons). Therefore, we put the following initial con- 
dition 

CEk(Ax) = C!k(Ax) = Cifl(Ax) = ^LRi^x) ^ afl(Ax) 

(6.6) 

Then, a model with Q{Ax) = <^5(Ax) = Q{Ax) is obvi- 
ously ruled out because we cannot give the observed val- 
ues of quark and charged lepton masses simultaneously. 
We must consider 



ej^(Ax) = ei(Ax)^a(Ax)^a(Ax) 



(6.7) 



We tentatively put ^|(Ax) = £,lr{Ax)- The mmrerical 
results are as follows: 

ai?.(Ax) = ^(Ax) = 1.20 , esi^x) = 0.80 , (6.8) 



As/Ar = 107 



(6.9) 
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zi = 0.01617 , Z2 = 0.2349 , = 0.9719 . (6.10) 

In the quark and charged lepton mass expressions (3.19) 
the factors and appear only in terms of the com- 
binations i^gAs and Cl^s, respectively, so that the abso- 
lute values of ^| and depend on the choice of the input 
value of As. Only the ratio Cs/Cs is substantial for the 
fitting of the quark and charg ed lepton mass. (However, 
as we state in the Sec. VIIl, the neutrino mass differ- 



ence between m^a and m,^2 rapidly varies in the range 
III. Therefore, in the neutrino mass matrix, the choice of 
the input value A5 is important.) We can obtain 



e|(A^)/C|(A^) 1.5 



(6.11) 



for any initial values of ^KAjf) with 0(1). The val- 
ues (6.10) are nearly in agreement with the values zi = 
0.01622, Z2 = 0.2357, and Z3 = 0.9717 which are obtained 
from Eq. (1.7) at ^ — mz- We can see that the effect of 
the evolution is not so large for Z^. 

The value of the parameter bci{Ax) is determined 
from the fitting of the observed down-quark mass ra- 
tios mj_/mg and rrig/mi, and the CKM matrix param- 
eter |Kis("^z)| ~ 0.22. In Fig. ||, we illustrate the mass 
ratios md{fi)/ms{fi) and ms{fi)/mi,{fj,) and the CKM 
parameter |Vus(/i)| at 11 — mz versus the parameters 
bd and (3d, where we have re-defined the complex pa- 
rameter bd by bde'^^'' with two real parameters. For 
convenience, in Fig. ^, the quantities are expressed in 
the unit of the corresponding observed values at = 
mz (for example, in Fig. ||, the curve md/ms denotes 
[md{fi)/ms{fi)]f,=mz/['md/ms]observed)- Wc obtaiu 



bd{Ax) - -1.20 , PdiAx) = 19.2° , 
which give the following predictions at — mz'- 



(6.12) 



mu{mz) 


= 2.60 X 10"^ 


GeV , 


mc{mz) 


= 6.92 X 10'^ 


GeV , 


mt{mz) 


= 182 GeV , 




md{mz) 


= 4.38 X 10^^ 


GeV , 


ms{mz) 


= 9.84 X 10^2 


GeV , 


mb{mz) 


= 3.02 GeV , 




me{mz) 


= 4.90 X 10"'' 


GeV , 


m^imz) 


= 1.03 X 10"^ 


GeV , 


mr{mz) 


= 1.76 GeV . 





(6.13) 



The experimental values corresponding to the results 
(6.13) are as follows pO[ |: 

m„(mz) = (2.331^5) x 10"' GeV , 

m,{mz) = (6.85t°;^?) x lO^^ GeV , 

mt{mz) = (181 ± 13) GeV , 

md{mz) = (4.691°:^°) x 10"' GcV , 

ms{mz) = (0.9341°;^^®) x 10"^ GeV , (6.14) 

mb{mz) = (3.00 ±0.11) GeV , 



meimz) = (4.8684727 ± 0.00000014) x 10"'' GeV , 
m^{mz) = (1.0275138 ±0.0000033) x 10"^ GeV , 
mr{mz) = (1.7467 ±0.0003) GcV . 

The results (6.13) is in agreement with the observed val- 
ues (6.14) within the experimental errors. 

The predicted values of \Vij \ depends on the phase pa- 
rameters 5( given by Eq. (1.4). Only when we take those 
as (1.8) (at /i = Ax), we can obtain reasonable values 
of \Vij\. For example, for 5^ — 5'^ = n, we obtain the 
predictions at fj. — mz 

iKsl- 0.220, lyehl = 0.0668 , 
\Vub/V,b\ - 0.0558 , 

\Vtd\ = 0.0177 , (6.15) 
J = 3.25x lO"'^ . 



The observed values |Tl[ are 

\Vus\ = 0.2196 ±0.0023 , 
\Vcb\ 0.0402 ±0.0019 , 
iKfc/Kfcl = 0.090 ±0.025 



(6.16) 



Although the results (6.15) are roughly consistent with 
experiments, the value \Vcb\ — 0.066 is somewhat large 
compared with the observed value \Vcb\ = 0.040. This 
discrepancy can be adjusted by considering a small devi- 
ation from TT of the relative phase S'!l; — S^ as demonstrated 
in Ref. ||. 

Related to the phcnomenological requirement (1.8), it 
is interesting to consider that which is the coefficient 
of the Higgs scalar (pL related to which is the coef- 
ficient of the scalar (j>L as 



(6.17) 



Then, the relations (1.8) mean that (i^/)ii and (y/)22 
are real, while (y/)33 is almost pure imaginary. We take 

(Z")t = Z'^ = diag(zi,Z2,^3e''^) . 



(6.18) 



The parameter 63 (= 6^ ~ ^S^) does not affect the 
masses, but only the CKM mixings. It is interesting to 
consider that the parameter Ss^Ax) takes its value such 
as the CKM mixings become minimum, i.e., such as the 
value \^iji-^x)\'^ takes the minimum. This require- 

ment gives the initial value (53(A3) = 84° (see Fig. ||). 
Then, we obtain the predictions of \Vij \ at /i = mz 



|Ks|= 0.220, l^efcl^ 0.0418 
\Vub/V,b\ - 0.0726 , 
\Vtd\ - 0.0109 , 
J = 2.38 X lO-'^ . 



(6.19) 



which is in excellent agreement with the experimental 
values (6.16). In Fig. H, we illustrate the predicted val- 
ues \Vij{mz)\ versus S^IAx)- As seen in Fig. |^, the value 
of 63{Ax) at which J2i^j takes the minimum 

also gives the minimum of the CKM mixings at /i = mz- 
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VII. NUMERICAL RESULTS IN THE SUSY 
MODEL 

The behavior of {^) in the SUSY model is somewhat 
different from that in the non-SUSY model. Since in the 
SUSY model, the top quark mass TOt(/u) is given by 



1 Vl 

^^^^^ " V^^^^^^V^^"^^^ ' 



(7.1) 



where vl/V2 = 174 GeV and tan/3 = tan/^L = Vl/""!: 
the initial value of Cti^^z) in the SUSY model corre- 
sponds to 



mmz) 



\SUSY 



= mmz)] 



1 



nonSUSY 



sin/3 



(7.2) 



However, this does not mean [£,i^{Ax)]susY — 
[^]^(Ajc)]„o„5[/5y/sin/3, because the behavior of 
[(,Lif^)]susY is considerably different from that of 
[iLit^)]no7iSUSY- In Fig. ||, we illustrate the behavior 
of f}^(^) in the SUSY model for the case of tan/3 — 3.5. 
If we take tan/3 < 2.5, the initial value of £,i,{mz) be- 
comes ^2(™z,tan/3 > 2.5) > ^^^(m^, tan/3 = 2.5) from 
Eq. (7.2), so that the curve of ^]^(At) will be illustrated 
in the upper side of the curve given in Fig. |. There- 
fore, for a case with a small value of tan /3, the Landau 
pole of ^£(/i) appears at a relatively lower energy scale. 
We consider that the model should be calculable per- 
turbatively, so that a case with such a large value of 
should be ruled out. As seen in Fig. ||, since the model 
gives, in general, CE(Ax) > CilA*) {mz < M < Ax), 
the value ll{Ax) should, at least, be [^l{Ax)?/4:n < 1, 
i.e., ^£(Ax) < ViT^ = 3.54. However, when we take 
contributions from the higher order corrections into con- 
sideration, even the value ^)^(Ax) — 3.0 is still dangerous. 
Therefore, we put the constraint ^]^(Ax) = 2.0 for the 
results of the present one loop calculation. In Fig. ^ we 
illustrate the predicted value of mt{mz) for the initial 
values S,liAx) = = 3.54 and CE(Ax) < 2.0, where 
we have used the input values 



Ax = 2 X IQis GeV , Ag 6 x 10^^ GeV 



(7.3) 



The value of A5 has been chosen as the neutrino mass- 
squared difference Amgj is of the order of (10^"^ — 10^^) 
eV2. 

From Fig. ||, we conclude that the value of tan /3 must 

be 



tan/3 > 3 



(7.4) 



Prior to the numerical investigation of the evolutions in 
the SUSY model, in order to see the difference between 
the parameter structures in the non-SUSY and SUSY 
models, let us give a rough sketch for the parameters in 
the case of the SUSY model by neglecting the evolution 
effects. The quark mass matrices M„ and Md are given 
by 



fc=i 



(7.5) 



fc=l 



(7.6) 



where = diag(l, 1), O'^ = diag(l, 1, 1/(1 -t- 3&d)), and 
Z is given by Eq. (5.14). Here, for simplicity, we have 
assumed [3l = — [3s = /3. For tan/3 > 3, the fac- 
tors sin /3 and cos /3 are approximated as sin /3 ~ 1 and 
cos /3 ~ 1/ tan /3, respectively. Obviously, the model with 
~ addition to the constraint 



e];fl(Ax)=^!«(Ax)^eL(Ax) , 



(7.7) 



is ruled out, because we cannot fit the up- and down- 
quark masses simultaneously due to the existence of the 
factor cos /3. Therefore, we must consider a model with 
7^ ?5 differently from the constraint (6.6) in the non- 
SUSY model. If we consider 



CI sin/?, el =^ el cos /3 



(7.8) 



then the model becomes similar to the case of the non- 
SUSY model, because 



A, 



sin/3 : 



— cos p , 



A, 



(7.9) 



and we will obtain reasonable fittings for the quark 
masses and CKM matrix parameters as well as in the 
non-SUSY model. Note that for a large value of tan /3, 
the value of K"^ = ^l^r/^s becomes large because K'^ ~ 
if" tan/? from the relation (7.9), so that we cannot eval- 
uate the RGE (3.12) perturbatively. We must take the 
value of tan/3 near to the lower bound given by Eq. (7.4). 

When we take the evolution effects into consid- 
eration, the situation is further complicated. The 
evolutions of ^{(m); ^^(m) and in the SUSY 

model are quite different from those in the non- 
SUSY model. We illustrate the behaviors of 
(fi) / (Ax) which correspond to the behaviors of 

K{(M)^^X/^)/e^(Ai)]/K{(A^)e^,(A^)/e^(Ax)] in the non^ 
SUSY model and those in the SUSY model in Figs. Q 
and H, respectively. In Fig. we see that the values 
TO„(/i) and mc(/i) cause rapid changes in the range HI. 
In the non-SUSY model, the charged lepton mass ratios 
are almost invariant, i.e., me(^)/m^(/x) ~ constant and 
TO^ (/i)/mT- (/x) — constant, while, in the SUSY model, 
the mass ratio (/i) /mr (/i) shows a considerable change 
(although me{^J.) ~ rUf^di) still holds). 

The situation is critical for the input values. If we 
adhere to the input value mt{mz) = 181 GeV, then it 
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is hard to obtain reasonable values of the other quark 
mass values mc, m„, rub, rris and md for any parameter 
values of As/Aji and bd- However, if we take a slightly 
lower value of mt{mz), for example, mt{mz) — 168 GeV 
[of. [mt{mz)]observed — 181 ± 13 GeV], we can find the 
following parameter values 

tan/3 = 3.5, As/A^ = 38 , (7.10) 
zi = 0.01449 , Z2 = 0.2117 , = 0.9772 , (7.11) 

Clr{^x) = CL(Ax) ^ eLRi^x) = 1.3 , ^£fl(A^) - 1.0 

(7.12) 

ej^(Ax) = 1.7, C|(Ax) = 0.50, C^(Ax) = 1.0, 

(7.13) 



are also invariant. Since the coefficients Hf^^ and H^j^ in 
the charged lepton sector are given by i?^^ — H'^j^ = 
in the ranges I and II for the non-SUSY and SUSY mod- 
els, the form of the charged lepton mass matrix Mg is also 
invariant below jj, — As- Therefore, the Maki-Nakagawa- 
Sakata (MNS) |jl^ matrix U — U^j^Ui, is invariant in the 
ranges I and II. Note that in the conventional model, the 
neutrino seesaw mass matrix can vary the from. The 
neutrino mass matrix in the present model can vary the 
form only in the range III (Ag < /i < Ax). The reason is 
that in the conventional model the scalar 0j couples to 
^L^Ri while that in the present model couples to Vj^En, 
so that the contribution of (pL to Hj^j^ in the latter case 
is decoupled below — As- 

For the numerical study, the case with bi, = ~l/2 
is most interesting, because the inverse matrix of 
Ysi^x) = ^S(Ax)[l + 3b4Ax)X] with b^iAx) = -1/2 
has the form 



bd = -1.2 , Pd = 19.4° 



(7.14) 



[Ys^i^x)Y 




^si^x) 



(8.1) 



which leads to the following quark and charged lepton 
masses and CKM matrix parameters: 



mu{mz) 


= 2.47 X 10"^ 


GeV , 


mc{mz) 


= 6.46 X 10"^ 


GeV , 


mt{mz) 


= 167 GeV , 




mdimz) 


= 4.49 X 10"^ 


GeV , 


ms{mz) 


= 1.00 X 10"^ 


GeV , 


mtimz) 


= 2.83 GeV , 




rrieimz) 


= 4.87 X 10"'' 


GeV , 


mf_,{mz) 


= 1.03 X 10"^ 


GeV , 


mr{mz) 


= 1.75 GeV , 





(7.15) 



|K.|= 0.220, |ye6|= 0.0665 
\Vub/V,b\ = 0.0603 , 
\Vtd\ = 0.0179 , 
J== 3.38 X 10"^ . 



(7.16) 



The values \Vij\'^ are again desirably adjustable by the 
phase parameter 6^ defined by (6.18). 



VIII. EVOLUTION OF THE NEUTRINO MASS 
MATRIX 

The evolution of the neutrino mass matrix — 
K" {(t)lf / {^^) is described by the RGB (3.13). Since 



the coefHcient H'^j^ in the ranges I and II is given by 
Xhl, (4.28), for the non-SUSY model, and by 
0, (A. 15) and (A.24), for the SUSY model, the 
form of the matrix if at ^ = A^ does not vary from that 
at II = As, so that the mass ratios and mixing matrix Ui, 



so that 



ZiZ2 Z2Z3 



The form (8.2) is well known as the Zee-type p^ ] 
matrix, which can lead to a large mixing [|l4| . 

The mass eigenvalues m^i and mixing matrix 
11 = Ax are given by 111 



m^3 



Z2Z3 - 



2Z2J 

—) 4 

2Z2 ' 



'0 I 



(8.2) 

mass 
U at 

(8.3) 



TO, 



1 



(a 



fAl 



/ 1 



U 



V 



- £1 

-zi 



V2 22 



l.£i(l 



1 



Z2) 



l.£l(l 



V2 ^2 



^2) 



V2 



1 
V2 



(8.4) 



(8.5) 



The model with bd = —1/2 gives highly degenerate mass- 
squared levels TO^2 — s-nd a large mixing between 
and Vr Sit fi = Ax- Therefore, the model has a possi- 
bility that it can give a reasonable explanation for the 
atmospheric neutrino data p^ . 

In Figs. ^ and ^ we illustrate the behaviors of the 



mass-squared differences Am 



in the non- 



SUSY and SUSY models, respectively. As seen in Figs. 
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and |l^, the mass-squared difference Arrigj rapidly in- 
crease according as the energy scale decreases in the 
range III. The numerical results are given in Table 0. 
We can see that the neutrino mass ratios are invariant in 
the ranges I and II. 

As we stated already, the values (zi, Z2, 23) (therefore, 
the mass ratios TOe/m^ and m^/mr) are almost invariant 
in the range III, while the ratio Am|2/Am|]^ is rapidly 
vary in the range III. Although the relations (8.3) give 
Am§2 — ^zlz2Zz{mQY and Ato2i ~ (22^:3)^(^0)^, the 
rapid decrease in the ratio Am32/ATO2i does not mean 
the rapid decrease in the ratio 2^/2223. The rapid de- 
crease comes from the slight deviation of the parameter 
hu{^J) from the value hv{^.x) = ^1/2. The value of 6/(/i) 
is not invariant in the range III, although the form of 
Yg (/i) , "the unit matrix plus a democratic matrix" , is 
invariant. When we denote 



bu{^i)^-^{l + e,{^l)) 

the expression (8.1) is replaced with 



1 /-2e. 1 1 \ 



(8.6) 



(8.7) 



so that 



-2e^zl 


Z1Z2 


Z1Z3 


Z1Z2 




Z2Z3 


ZlZs, 


Z2Z3 


2s,yZ^ 



Therefore, the mass eigenvalues in the range III are given 

by 



m^i ~ -2(1 + 3£y)z^TOo , 



m^3 



Z2Z3 
Z2Z3 + 



(8.9) 



2^2 / 



instead of (8.3), and the mass squared differences Am^i 
and Arn32 are given by 

Am2i ~ (22^3)^(^0)^ , 



Am32 ~ Az2Z3{zl ~ e^)(mo)^ 



(8.10) 



Note that the approximate expression (8.19) tell us that 
Am|2(A*) takes a zero between /i = Ax and /i = A5 be- 
cause ey{kx) = < Zi(Ax) z\{ks) < e^{As), e.g., 
e^As) = 7.3 X 10-2 and z^{Ax) =i zf{As) ~ 2.6 x lO"'' 
for the non-SUSY and e^{As) = 1.1 x 10"^, zf{Ax) =i 
2.1 X 10-4 and zf{As) ~ 2.6 x 10"'' for the SUSY model. 
In fact, we can see this at a point which is very close to 
/X = Ax in Figs, ^and ^ Thus, the value of Ar7i32(^) is 
highly sensitive to the value of £,y(^), although Ato2i(m) 
is not so. 



In general, since the mixing angle ^23 is given by 

2(M,)23 



sin 26*23 



m^3 - m^2 



(8.11) 



the mixing angle ^23 in the conventional democratic type 
neutrino mass matrix model is sensitive to the energy 
scale , because Am§2 (a^) has a large energy scale de- 
pendency. In contrast to the conventional model, the 
mixing angle ^23 in the present model does not so dras- 
tically vary. The reason is as follows: the neutrino mass 
matrix Mi, in the present "democratic" seesaw model is 
not democratic, i.e., the form of is given by Eq. (8.2). 
In fact, the present model gives not 771,^2 — ^1,^3, but 
m^2 — — w,y3, so that the evolution effect on J723 is not 
so sensitive as seen in Eq. (8.11). 

As seen in Table ^ the model can fit the value Am§2 
to the atmospheric neutrino data ||l6| {A'm'^2) observ = 
3.2 x 10-'^ eV^ by adjusting the value of A5, but it can- 
not give any explanation of the solar neutrino data , 
because of Am|;^ ^ Am|2 = Am^j^. We must intro- 
duce a further mechanism for the explanation of the so- 
lar neutrino data, for example, as discussed in Ref. 
However, since the purpose of the present model is not 
to propose a plausible neutrino mass matrix model in the 
framework of the USM, but to see the characteristic fea- 
tures of the neutrino mass matrix evolution in contrast 
to the conventional seesaw model. Therefore, we do not 
touch the numerical fitting furthermore. 



IX. CONCLUSIONS 

In conclusion, we have investigated the evolutions of 
the quark and lepton mass matrices Mf (/ = u, d, v, e) 
in the universal seesaw model with detMp = in the 
up-quark sector F = U. 

The assumptions which have made in the present paper 
are classified into the following three categories: 
(A) Basic assumptions in the universal seesaw model with 
detM,7 = 0; 

B) Basic assumptions in the democratic seesaw model 
PI (we have taken the model as a more concrete one of 
the universal seesaw model with detA'% = in order to 
give an explicit evaluation of the universal seesaw model); 
(C) Tentative assumptions for convenience of the numer- 
ical evaluation. 

The assumptions in the category (A) are as follows: 
(Al) SU(3)^ X SU(2)^ X SU(2)^ x U(l)y. x \]{l)x gauge 
symmetries with the symmetry breaking pattern (2.1); 
(A2) Hypothetical heavy fermions F = {U,D,N,E) 
which belong to (1,1) of SU(2)^ x SU(2)^ and acquire 
masses of the order of A5 at the energy scale ^ = A5 
except for U3L and U3l- 

In the present model, therefore, the quark and charged 
lepton mass matrices Mf (/ = m, d, e) and neutrino mass 
matrix are given by 
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Mf 



M,^Yl{Yl^)-\Ylf{Kl/Ks) , 



(9.1) 



(9.2) 



except for the top quark, where = ^-R = 

and As = ($°). The evolutions below = A5 are de- 
scribed by the RGE (3.12) and (3.13) for the seesaw op- 
erators. On the other hand, the top quark mass mt{fi) 
given by the expression (3.18) is still described by RGE 
(3.6) for the Yukawa coupling constants below /i — A5. 
Although the heavy fermions F do not contribute to the 
evolutions below /x = A5, the third family "would-be" 
heavy up-quark U3 can contribute to the RGE even be- 
low n = As- However, as far as the i?^^ {F — v^e) 
and iJ^fl. terms in the lepton sectors are concerned, the 
would-be heavy quark cannot contribute to those, so 
that the forms of the mass matrices M^{^) and Me{ii) 
are invariant below fi = A5. 

The assumptions in the category (B) are as follows: 
(Bl) At a unification scale /i = Ax, the Yukawa cou- 
pling constants f/ and have the same form, i.e., 

y/(Ax) = ri(Ax)^r/^(Ax). 

(B2) At II — Ax, the heavy fermion mass matrices Mp 
(therefore, the Yukawa coupling constants Yg) [and also 
the Majorana masses Ml {Mr) of the neutral fermions 
Nl {Nr)] take a simple diagonal form (5.1), the form " 
the unit matrix pulse a democratic matrix" , on the basis 
on which the Yukawa coupling constants Y^^(Ax) are 
diagonal. Then, the form (5.1) is invariant under the 
evolution in the range III. 

(B3) The values of the parameter &/ in the matrix Yg 
given by Eq. (5.1) are given by hf, = 0, = —1/2 and 
bu — —1/3 a.t fi — Ax- (The value fe^ is kept as a free 
parameter in order to fit the up- and down-quark masses 
and CKM matrix parameters reasonably.) 

In this model, the top quark mass mt(/i) is given by 
(5.15). The behavior of mi(/i), i.e., f}^(/i), is given in 
Figs. |l| and | for the non-SUSY and SUSY models, re- 
spectively. We can obtain the constraint on the values of 
the intermediate energy scales Ar and As by considering 
that the model should be calculable perturbatively. In 
the non-SUSY model, since As/Ar ^ 10^ from the ratio 
mt/rric, we find the constraint 



10^" GeV < As < 10'^ GeV 
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(9.3) 



for Ax 10^^ GeV. In the SUSY model, the results 
highly depend on the input parameter tan /3. From the 
numerical study, we obtain the constraints 



3 < tan /3 < 4 , 



10^° GeV < As < 10'^ GeV , 
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(9.4) 



(9.5) 



for Ax ~ 10^^ GeV. (The above numerical results are 
slightly dependent on the assumptions stated below, but 
the dependence is not so large.) 



The assumptions in the category (C) are as follows: 
(CI) For convenience of the numerical evaluation, 
52l(/^) — 52_r(a') has been assumed. Then, we can as- 
sert Y[{p) = Y^{ii) in the range III (As < ^ < Ax) as 
we have shown in Eq. (4.8). 

(C2) For evaluation of the non-SUSY model, the initial 
condition 

Clr{^x) = e!«(Ax) = ei«(Ax) = eL«(Ax) (9.6) 

has been assumed together with the initial condition 
(3.5), i.e.. 



(9.7) 



However, since there is no solution of the parameter val- 
ues for the SUSY model under such a constraint (9.6), 
the constraint corresponding to (9.6) in the SUSY model 
has been loosened as 

C£«(Ax) = C!«(Ax) ^ eii^(Ax) - ilR{Ax) , (9.8) 

although the initial condition (9.7) has still been required 
in the SUSY model. 

(C3) For the non-SUSY model, we have assumed 
?s(Ax) = i%{Ax) ^ ei(Ax) = eS(Ax), but, for the 
SUSY model, we have assumed that each value of (,^g{^ 
may be different among them, because the previous 
condition is too strong for the SUSY model and the 
up<-!-down symmetry is already broken due to the fac- 
tor tan /3 ^ 1 in the SUSY model. 

In the conventional model for quark and charged lep- 
ton masses (i.e., not seesaw model), the following approx- 
imate relations are satisfied in the non-SUSY and SUSY 
models: 

{mulmc)L {md/ms)L {rne/m^)L {m^_,/rnr)L 



(m„/TOc)x {■md/ms)x {me/m^j)x {m^_,/mr)x 

(9.9) 



{mu/mt)L {mc/mt)i 



{mu/mt)x imc/mt)x 



{md/mbjL (ms/mb)! 



{md/rrib) 



X 



{ms/mb)x 



1 + Ed , 



(9.10) 



(9.11) 



\Vcb{AL)\ |Kb(AL)| \Vtd{AL)\ 



\Vcb{Ax)\ |K6(Ax)| \Vtd{Ax)\ 



l + Sd, (9.12) 



where (rau/mc)L denotes m„(AL)/mc(AL), and so on. 
The relations (9.9)-(9.12) are due to that the Yukawa 
coupling constant yt of the top quark in the conventional 
model is very large compared with the other Yukawa 
coupling constants. In the present model, as seen in 
Figs. and |, the relations (9.9)-(9.12) are also satis- 
fied in the range I (Al < fJ. < Ar) (so that we read the 
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relations (9.9)-(9.12) as X — !■ R). The values of e„ and 
Ed are approximately given by e„ ~ for the non-SUSY 
model, and by ~ Ssd for the SUSY model. In the 
range II (A^^ < /x < As), the relations (9.9)-(9.12) are 
slightly broken. In the SUSY model, the values show not 
Eu — Ssd, but e„ ~ Ed in the range II. However, in the 
model with Al/Ah ^ 1, which is required in order to 
make the neutrino masses tiny, the evolution effects in 
the range II are not so large, so that we can regard that 
the relations (9.9)-(9.12) are still satisfied in the range 
Al < < As, i.e., 

i?„(Ai) ~ ^^4t4(1 + - euS)DUAs) , (9.13) 



Dd{AL) 



TOfc(As) 



il + Ed)il-EdS)Dd{As) , (9.14) 



VIAl) ~ (1 + EvS)V{As)il + EvS) - 2EvVtdS , (9.15) 

where Du = diag(m„, mc, mj), Dd = diag(mrf, m^, mb), 
and S is defined by Eq. (4.9). In the present model, the 
value of Ey is not always given by ev — Ed because of the 
presence of the range II. 

Also in the ranges I and II, differently from the con- 
ventional seesaw model (for example, see Ref. |^), the 
neutrino mass ratios and mixing angles are not affected 
by the evolution effects: 



is not sensitive to the energy scale. In contrast to the 
conventional seesaw neutrino mass matrix, note that the 
present neutrino mass matrix M^, is form-invariant below 
fi ~ As, so that the neutrino mass ratios and mixings are 
invariant below /i — Ag. 

In the present paper, we have assumed SU(3)^ x 
SU(2)^ X SU(2)^ X U(l)rg X U(l)^ symmetries above 
/i = As. As seen in Figs.^and ^ in general, the rapid 
increasing of the Yukawa coupling constant Y^{iJ,) causes 
above n = As, although we have been able to find a set of 
the reasonable parameter values without having the Lan- 
dau pole below /i = Ax- The rapid increasing is mainly 
due to the rapid increasing of the gauge coupling con- 
stant gi above /i — As- If we want to build a unification 
model with a unified gauge symmetry G, we may consider 
that the U(l) symmetry is embedded into the unified 
symmetry G. (For example, see an SO(IO)^ x SO(IO)^ 
model |20), where SO(IO)^ x SO(10)j^ is broken into 
[SU(2) X SU(2)' X SU(4)]l x [SU(2) x SU(2)' x SU(4)]fl.) 
Then, the gauge structure above fi = As is different from 
the present model, so that the evolutions will be also 
different from the present results. (Of course, the evo- 
lutions below /i = As are still the same as those in the 
present paper.) It is likely that the gauge structure above 
fi = As is different from the present model. Our next 
task is to investigate what gauge structure above fi = As 
is promising for a unified description of the quark and 
lepton masses and mixings. 



mui{AL)/m^j{AL) 
m^i{As)/mi,j{As) 



V^j-(Al) 
Vij{As) 



1 



(9.16) 



(9.17) 



Note that the relation (9.16) does not mean 
Am|^(Ai)/Am|^(As) ~ 1. However, the ratio 
Am2x/Am§2 is again invariant in the ranges I and II. 

In the range III (As < fi < Ax), the relations (9.9)- 
(9.12) [(9.13)-(9.15)] and (9.16)-(9.17) are not satisfied 
at all. For example, the behavior of A?7i§2(m) is highly 
sensitive to the value £i/(/x) and is given by Eq. (8.10). In 
other words, the differences of the numerical behaviors of 
the quark masses, CKM matrix parameters and neutrino 
mass squared differences from those in the conventional 
model are substantially formed in the range HI. 

Note that the mass ratios me/nifj, and mu/nic are 
almost constant (although the ratio niu/mc is slightly 
changed in the SUSY model), so that the phenomenolog- 
ically well-satisfied relation (1.14) still holds under the 
evolutions. 

For the neutrino mass matrix M^,, we have investigate 
the model with &^(Ax) = —1/2, which leads to a large 
mixing sin^ 023 — 1- Although the mass-squared differ- 
ence Am32(/i) is highly sensitive to the energy scale /x 
in the range III (As < < Ax), the mixing angle 623 
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Appendix 

In Sees. 4 and 5, the coefficients of RGE (3.6), (3.12) 
and (3.13) have been given only for the case of non-SUSY 
scenario with one SU(2)-doublet Higgs scalar. In the 
present Appendix, we give the coefficients of RGE in the 
minimal SUSY scenario. 
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where A = L.R. 
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where A — L, R and f — u,d, e. 
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TABLE I. Quantum numbers of the fermions / and F and Higgs scalars cpL, 4'R and $ for 
SU(2)^ X SU(2)^ X U(l)^^ X U(l)^. 
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TABLE IL The squared mass difference Amfj = mjj — m^^. The values of the input parameters are the same as in Figs. 9 
and 10. The absolute values of |Am|j| should not be taken rigidly, because we can adjust those by the value of As. 





non-SUSY model 


SUSY model 




at jj, = Al 
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at /J = Ax 
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at n = Ax 


\AmU [eV^] 
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4.08 X 10"" 


lAm^il [eV^] 
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1.01 X 10"^ 


lAmi^a/Am^il 
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0.00484 
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FIG. 1. Behaviors of CM in a non-SUSY model for 
the cases (a) As = 10® GeV, (b) As = 10« GeV, (c) FIG. 3. J2i^j\Vij{^x)\^ versus 5i{Ax). 

As = 10^^ GeV, and (d) As = 10^^ GeV. The input values 
are mt{mz) = 181 GeV and As/Ar = 107. 




FIG. 2. Predictions of md/ms, ms/mb and |Vii.s|, and 
their dependency on the parameters hd and (id- Here, the 
mass ratios are denoted in the unit of the corresponding ob- 
served values which are quoted from Ref. [10]. The dashed, 
solid and dotted lines denote ba = —1.1, —1.2 and —1.3, re- 
spectively. 



FIG. 4. Predicted values of the CKM matrix parameters 
|Vij(mz)| versus the parameter S^{Ax)- Other input values 
of the parameters are As = 3 X 10^^ GeV, As/Ar = 107, 
bd = -1.2 and Pd = 19-2°. 
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FIG. 5. Behaviors of ^"(p) in a SUSY model for the cases 

(a) As = 10® GeV, (b) As = 10^ GcV, (c) As = lO^^ GeV, 
and (d) As = 10^^ GeV. The input values are mt(mz) = 181 
GeV, As/Ail = 38 and tan/3 = 3.5. 



FIG. 7. Behaviors of mj {fi)/m{ (Ax) (/ = u,d,v,e; 
i = 1,2, 3) in the non-SUSY model. The dotted, broken and 
solid lines denote the first, second and third fermion masses, 
respectively. The input parameter values are As = 3 x 10^^ 
GeV, As/Ah = 107 and 6d(Ax) = -1.2e"^-2°. 
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FIG. 6. The top-quark mass mt{mz) versus tan/3 in 
a SUSY model. The solid and broken lines denote the 
cases with the initial conditions (a) Cl(Ax) ~ 2.0 and (b) 
i^i^x) = \/47r = 3.54, respectively. The other input values 
are As = 6 x lO" GeV and As/Ar = 38. The horizontal 
solid and broken lines denote the center and lower values of 
the observed top quark mass aX n — rnz, respectively. 
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FIG. 8. Behaviors of ml [fx) / ml {Ax) (/ = u,d,i',e; 
i = 1,2,3) in the SUSY model. The dotted, broken and 
solid lines denote the first, second and third fermion masses, 
respectively. The input parameter values are As = 6 x 10^^ 
GeV, As/Ar = 38, tan/3 = 3.5 and 6ci(Ax) = -1.2e"^-*°. 
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FIG. 9. Behavior of |Am-j(/x)| in the non-SUSY model. 
The input parameter values are the same as in Fig. 7 with 

a = a {A = L,R,s). 



FIG. 10. Behavior of \Amij{n)\ in the SUSY model. The 
input parameter values are the same as in Fig. 8 with = 
{A = L,R, S). 



21 



